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ABSTRACT 



The spectra and decay rates of cc and bb levels are well described, for the most part, 
by a power-law potential of the form V(r) = X(r a — l)/a + const., where a ~ 0. 
The results of an up-to-date fit to the data on spin-averaged levels are presented. 
Results on electric dipole transitions in systems bound by power law potentials are 
also presented, with applications to the bottomonium system. 

Charmonium (cc) and bottomonium (bb) systems provide a rich source of 
information on the interquark force at distances ranging from less than 0.1 fm to 
greater that 1 fm. At short distances our theoretical prejudices favor a potential 
which should act like a Coulomb potential V(r) = a s (r)/r, with a s becoming smaller at 
shorter distances owing to the asymptotic freedom of QCD. 1 At long distances, there 
are both experimental and theoretical reasons 2 to believe that the interquark force 
in QCD becomes approximately distance-independent, corresponding to a linear 
potential V(r) ~ r. The cc and bb systems appear to lie in an intermediate range, 
where a power-law potential V(r) — r a provides a convenient interpolating form 3 
between the short-distance Coulomb-like and long-distance linear behavior. 

An early fit to quarkonia spectra 4 found a power a ~ 0.1. Since then, data on 
P-wave levels have appeared, 5 6 and information on leptonic widths has improved. 
It is appropriate to update the earlier fit for a number of reasons. The power-law 
method can be of use in estimating properties of systems containing b and c quarks, 
and in interpolating between these cases to get estimates of be properties. The power 
law method also gives an estimate of the mass difference between the b and c quarks, 
which can be of use in attempts 7 to extract the Cabibbo-Kobayashi-Maskawa matrix 
element V C b from data on semi-leptonic b decays. One would also like to see if there 
is a consistent pattern of data signaling a departure from a single effective power 
at short distances, as one would expect from the short-distance behavior of QCD. 
Finally, power law potentials can be of use in explaining patterns of electric dipole 
transition rates in the bb system. 

Since the power-law description does not give an adequate description of 
spin-dependent effects, we fit spin-averaged levels. For S waves, we use 

M(S) = [M^ So) + 3Af( 3 5!)]/4 , (1) 
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a combination which eliminates hyperfine splittings. For P waves, we use 

M(P)=M( 3 P) = [M( 3 P ) + 3M( 3 P 1 ) + 5M( 3 P 2 )}/9 , (2) 

which eliminates spin-orbit and tensor force splittings. In the latter case we assume 
hyperfine effects are small. In the case of bb levels, the masses of the X S levels are 
not experimentally known. However, the hyperfine splitting can still be estimated 3 
using information from the leptonic widths. 

We obtain theoretical values for the masses and leptonic widths of the levels 
using non-relativistic quantum mechanics. We find the energies and wavefunctions 
of the radial Schrodinger equation by solving the dimensionless equation numerically 
and then rescaling the dimensionless quantities by appropriate powers of the mass 
and coupling constant. 3 We take V(r) = \(r a - l)/a + C. Particle masses are then 
given by M = E + 2toq, where E is the binding energy and tuq is the quark mass. 

We present results from three fits. In each case we minimize x 2 , with the 
standard deviations for the energies set equal to 10 MeV, and the standard devia- 
tions for the leptonic widths set equal to the experimental errors. In the first fit, we 
consider the levels only. In this case, we find that the best fit is given by a potential 
V(r) ~ r -o o45^ anc j that the quark mass difference is m b - m c = 3.19 GeV. However, 
the fit gives little preference for quark masses, and in fact the best fit is given by 
m b , m c — > oo. In the second fit, summarized in Table I, we remedy this by including 
the leptonic widths, which are given by the formula 8 

_ 16a s {m Q ) 

and where we use 9 a s (m b ) = 0.189 ±0.008, a s (m c ) = 0.29 ±0.02. In this case, we find 

m b = 5.24 GeV , m c = 1.86 GeV , m b - m c = 3.38 GeV (4) 

for the quark masses, while in the potential V(r) = \{r a — l)/a+ C, we have 

a = -0.14 , A = 0.808 , C = -1.305 GeV . (5) 

We find that the fitted masses are correct to within 10 MeV or so, but the leptonic 
widths are off, typically, by 20 to 30 percent. In an attempt to remedy this, in the 
third fit, we include an ad hoc relativistic correction to the leptonic widths. We 
correct Eq. (||) by introducing a factor of the form (1 + K(v 2 /c 2 )), and treat K as 
a free parameter. In this case, we find a = -0.12, A = 0.801, and C = -0.772 GeV, 
while for the quark masses we have m b = 4.96 GeV, m c = 1.56 GeV, and m b - m c = 3.40 
GeV. The constant K is found to be 1.25. The fitted quantities are given in Table 
I. The relativistic correction gives only a marginal improvement in the fit. In the 
fits to masses and leptonic widths, we find that m b and m c are not particularly well 
determined. However, the mass difference m b - m c appears to be more stable, with 
a value around 3.39 GeV. 

We can use the fitted potentials to estimate the centers of gravity of a few 
low-lying levels of be, bs and cs levels. For the be system, we find 



r(QQ^e+e-)= ^ j*(0)|' 



(3) 



M b5 (lS) = 6.304 GeV, M bs (2S) = 6.898 GeV, M be {lP) = 6.764 GeV, 



(6) 



Table 1: J/ip and T masses and leptonic widths 



Particle 




Mass (GeV) 






Width (keV) 






(Expt.) 


(NR a ) 


(RC b ) 


(Expt.) 


(NR a ) 


(RC b ) 


J/V(1S) 


3.068 


3.077 


3.079 


5.36±0.29 


6.41±0.43 


6.29±0.43 


J/V(2S) 


3.663 


3.654 


3.654 


2.14±0.21 


2.03±0.13 


2.04±0.13 


J/V(1P) 


3.525 


3.524 


3.522 


— 




— 


T(1S) 


9.449 


9.420 


9.423 


1.34±0.04 


1.21±0.02 


1.18±0.02 


T(2S) 


10.018 


10.044 


10.042 


0.563±0.14 


0.477±0.009 


0.475±0.009 


T(3S) 


10.351 


10.358 


10.358 


0.44±0.07 


0.285±0.006 


0.284±0.005 


T(4S) 


10.578 


10.564 


10.567 


0.24±0.05 


0.197±0.004 


0.200±0.004 


T(1P) 


9.900 


9.903 


9.900 








T(2P) 


10.260 


10.269 


10.267 








T(1D) 




10.181 


10.177 








T(2D) 




10.436 


10.435 









a No rclativistic corrections "With relativistic corrections 



while for the bs and cs levels, we have 

M CS (1S) = 2.085 GeV, M cS {lP) = 2.509 GeV, M bs (lS) = 5.401 GeV. (7) 

The experimental spin-averaged masses 10 ' 11 of the cs states are 2.075 GeV for the IS 
level, and 2.536 GeV for the IP level. We have also estimated the IS - 2S splitting 
in toponium. Taking m t =130 GeV, we find the splitting to be roughly 0.8 GeV. We 
expect that this is a conservative lower bound, since the short-distance behavior of 
QCD (which will make this splitting larger) should be important in this case. 

Power law potentials also offer some insight into the strengths of various El 
transitions in the bb system. 12 Experimentally, 10 we find that the 3S-1P transition 
is suppressed relative to the 3S-2P, and that the 2P-1S transition is suppressed 
relative to the 2P-2S. We can gain some insight into this by considering the the 
radial dipole matrix elements (unt\r\u n 'e±i), where n denotes the number of nodes in 
the radial wavefunction u n i. We can approximately evaluate these matrix elements 
by considering the large I limit. In this case, we expand the potential in the radial 
Schrodinger equation about the point the point f = £ 2 /( 2+Q ). The potential then has 
the form of a harmonic oscillator potential, plus anharmonic terms which we include 
perturbatively. For the dipole matrix elements, we find 

{u n z\r\uni-i) = £ 2 ^ a+2 \ (tW-ilrK-if) = ^*+(aK (2 ~ Q)/(4+2Q \ 

(un-ii-xMiw) - V^MaK (2 ~ Q)/(4+2a) , (uo*|r|ua<_i) = $(a)^ a/(2+a) , (8) 

where *± and $ are functions which depend only on a. We see from Eqs. (|S|) that 
transitions with An = are dominant, and that the others are suppressed by a factor 
£-A n /2 j n limit of large I. Furthermore, in the region of interest for quarkonium, 
a ~ 0, the coefficient $ is quite small: $(0) ~ -0.04. This further suppresses the El 
rate for An = 2 transitions. Direct numerical calculations show that this is also the 



case for small t we find that the 3S-1P rate is dramatically suppressed relative 
to the 3S-2P rate. Consequently, it appears that the suppression of El transitions 
with nonzero An is a general property of power law potentials. 

We can use the ratios of radial dipole matrix elements for certain transitions 
in the bb system to place limits on the power a. From the experimental data 610 we 
deduce the ratios 

Computing these ratios numerically and comparing with experiment, we find that a 
is constrained to lie in the region -0.2 < a < 0, consistent with what was found from 
the fit to levels and leptonic widths. Finally, we can estimate m b using the rates 
of the 3S-1P, 3S-2P, and 2S-1P electric dipole transitions. Extracting the radial 
dipole matrix elements from the data, we find that the potential model favors a 
value m b ~ 4.0 ± 0.9 GeV. This is smaller than expected, and in particular we find 
that the 2S-1P transition yields the surprisingly small value m b = 3.45 ± 0.44 GeV. 
This is an indication that either the potential model is inadequate for describing 
El rates, or that the T(2S) total width has been overestimated. 
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